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Abstract 

We consider long-range interactions between two alkali-metal atoms in their respective ground 
states. We extend the previous relativistic many-body calculations of Cq dispersion coefficients 
[Phys. Rev. Lett. 82, 3589 (1999)] to higher-multipole coefficients Cs and C\q. A special attention 
is paid to usually omitted contribution of core-excited states. We calculate this contribution within 
relativistic random-phase approximation and demonstrate that for heavy atoms core excitations 
contribute as much as 10% to the dispersion coefficients. We tabulate results for both homonuclear 
and heteronuclear dimers and estimate theoretical uncertainties. The estimated uncertainties for 
C§ coefficients range from 0.5% for Li2 to 4% for Cs2- 
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I. INTRODUCTION 



We carry out accurate relativistic many-body atomic-structure calculations of van der 
Waals interactions between alkali-metal atoms in their respective ground states. These 
long-range interactions may be parameterized using dispersion (van der Waals) coefficients 

C n 

V (R)a-Q-°*-°™+... (1) 

where R is the internuclear separation. A renewed interest in high-accuracy interatomic 
potentials has been stimulated by advances in studies of ultracold collisions 2]. At low 
energies, collision properties are typically very sensitive to details of the potentials. Thus 
accurate potentials are essential for reliable ab initio description of ultracold collision prop- 
erties and, conversely, a wealth of information about the potentials may be inferred from 
photoassociation and Feshbach-resonance spectroscopy with ultracold atomic samples. In 
particular, only recently interpretation of experiments with ultracold atoms allowed several 
groups to reduce uncertainties in the Cg coefficients to a fraction of a per cent 0,111, These 
inferred coefficients are in an excellent agreement with our values predicted using many-body 
perturbation theory jfjj. Even more refined understanding of details of ultracold collisions 
led very recently to constraints on higher-multipole coefficient C 8 for Rb_0, jl and Cs EJ. 
This latest progress and discrepancies between previous determinations [loT fill H2I . [lj| 
of C 8 and C w coefficients motivate us to calculate these coefficients using accurate relativis- 
tic many-body techniques of atomic structure. In particular, we demonstrate that usually 
omitted contribution of core-excited states increases C n for heavy atoms by as much as 10%. 

The main result of the paper — compilation of van der Waals coefficients Cs and Cio for 
homonuclear and heteronuclear Li, Na, K, Rb, and Cs dimers is presented in Tables IVIl 
IIXI The rest of the paper is organized as follows. In Section |H] we present the formalism. 
Numerical evaluation is discussed in Section ILTT1 A detailed analysis of underlying multipole 
dynamic and static polarizabilities is presented in Section IIVI Finally, in Section IS we 
compile dispersion coefficients and estimate theoretical uncertainties. Atomic units (|e| = 
m e = h = 1) are used throughout the paper. 



II. GENERAL FORMALISM 

The long-range part of electrostatic interaction between two atoms a and b in their re- 
spective spherically-symmetric states may be represented as 

V{R) = -Y J Ct/R 2 \ (2) 

n=3 

where R is the distance between atoms. For ground-state atoms van der Waals coefficients 
are given by j3| 

ct = £ I of (M «f M*- • (3) 

where I' = n — I — 1; af(iu) and a b v {iuj) are, respectively, 2 z -pole dynamic polarizability 
of atom a and 2 l -pole dynamic polarizability of atom b. The dynamic polarizabilities in 
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Eq. are defined as 



ai{iuj) = 2 Re 2^ j-i \ , ■ • (4) 



k 



(E k - E ) + iu 



Here the summation extends over a complete set of atomic states and are the zeroth 
components of spherical tensors of electric-multipole operators 



N 



ri^E^cffC'O, (5) 



=i 



where are normalized spherical harmonics |l6( and the sum is over all N atomic electrons. 

Previously many-body calculations of dispersion coefficients Cq were carried out in 
Refs. 0, [l?]], and here we focus on dispersion coefficients C 8 and C w . As follows from 
an examination of Eq. (J3J), we need to compute dipole ai, quadrupole «2, and octupole «3 
dynamic polarizabilities. In this work we employ dynamic dipole polarizabilities calculated 
previously in Ref. and determine higher-multipole polarizabilities «2 and a-^. 

Following we separate all intermediate states in the sum Eq. (@J) into valence and 
core-excited states 

cti(iu) = ajftiw) + ot1{iuj) + af°(iu) , (6) 

Here afiiu) is a traditional term encapsulating excitations of the valence electron. Con- 
tributions of electric-multipole excitations of core electrons are denoted by af(iu). Finally, 
a small counter term af{iuj) is related to excitations of core electrons to occupied valence 
state. We include these exclusion-principle-forbidden excitations in the calculations of core 
polarizabilities and thus we have to introduce the counter term (see Ref. [HI for more de- 
tails). We will discuss calculations of the ctf{iuS) and af(iu) terms later on. Here we just 
briefly comment on the counter term af{iuj). For octupole polarizabilities oef(iu) term sim- 
ply vanishes in independent-particle approximation since E3 selection rules would require an 
excitation from / shell to valence s-state and none of the alkalis considered here (Li through 
Cs) has filled f-shells. Since we employ dipole polarizabilities from Ref. 0, the counter 
term, calculated in Dirac-Hartree-Fock (DHF) approximation is included in a\{iu). Finally 
we disregard this correction for quadrupole polarizabilities, it gives a negligible contribution 
due to required excitation of deeply bound d electrons from the core. 

High-accuracy calculations of the dipole dynamic polarizabilities were carried out ear- 
lier in Ref. 0] and we employ these dipole polarizabilities in the present work. In those 
calculations a combination of several relativistic many-body techniques was employed. A 
dominant contribution to a\ has been calculated with all-order linearized coupled-cluster 
method truncated at single and double excitations. High-accuracy experimental values for 
energies and electric-dipole matrix elements for principle transitions has been employed to 
refine the dipole polarizabilities. In the following we focus on the quadrupole and octupole 
polarizabilities. 

To find the quadrupole a" 2 and octupole valence contributions we applied a relativistic 
many-body method initially suggested in Refs. jl^, H3| and subsequently developed in j^H 
221 . 123L Eil j. In this method one determines wave functions from solution of the effective 
many-body Shrodinger equation 

H&(E n )\V n )=E n \V n ), (7) 
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with the effective Hamiltonian defined as 



H cS (E) = H FC + £(£) . (8) 

Here Hpc is the frozen-core Dirac-Hartree-Fock Hamiltonian and self-energy operator £ 
is the energy-dependent correction, involving core excitations. Qualitatively E operator 
corresponds to core polarization term in model potentials employed in Refs. 0,13. In 
the present calculation the self-energy operator recovers second order of perturbation theory 
in residual Coulomb interaction and additionally accounts for certain classes of many-body 
diagrams in all orders of perturbation theory. 

The concept of effective Hamiltonian H e ^ may be extended to other operators. We 
introduce effective (or dressed) electric-multipole operators T l eS acting in the model space 
of valence electrons. These operators were obtained within the relativistic random-phase 
approximation (RRPA) jUl HH, E3] • Qualitatively, the RRPA describes a shielding of the 
externally applied electric-multipole field by the core electrons. The RRPA sequence of 
diagrams was summed to all orders of the perturbation theory. 

Once the ground-state wavefunctions are obtained from Eq. (JJJ), the dynamic valence po- 
larizabilities a\(ioS) are computed with the Sternheimer ji?} or Dalgarno-Lewis |28[ method 
implemented in the DHF+E+RRPA framework. (In the following we denote E+RRPA cor- 
rections as the many-body perturbation theory (MBPT) corrections.) Given ground-state 
wave-function \l/o and energy Eq, we find an intermediate-state wave function from an 
inhomogeneous equation 

l*/> = Re (- 1 —— £l^><M r o)eff|*o}} 

-"off — Eq + i J 

h VT^^^}- (9) 

-n e ff — E +IU) ) 

With such introduced ^ff Eq. (J3J) becomes simply 

aftiu) = 2{* \(I*U\* f ), (10) 

where subscript v emphasized that only excitations of the valence electron to higher virtual 
orbitals are included in the intermediate-state wave function \1/ f due to a presence of H e g in 
Eq. (jnj). As to additional contribution a\ of core-excited states, we employ the relativistic 
random-phase approximation method described in Refs. ji^, |26| . 



III. DETAILS OF NUMERICAL CALCULATION 

At the first stage of calculations we determined core orbitals and valence orbitals for 
several low- lying states from the frozen-core Dirac-Hartree-Fock equations |2fll |. The virtual 
orbitals were determined with the help of a recurrent procedure |30j . One-electron basis sets 
of the following sizes were used on the stage DHF+E calculations: 

Li : 1 - 17s, 2 - Yip, 3 - 16d, 4 - 16/, 5 - 10#; 

Na : 1 - 18s, 2 - 18p, 3 - lid, 4 - 17/, 5 - llg; 
K : 1 - 19s, 2 - 19p, 3 - 18d, 4 - 19/, 5 - \2g; 
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TABLE I: Comparison of DHF and many-body one-electron removal energies E va \ for Cs with 
experimental values. E va \ are given in atomic units. A are excitation energies from the ground 
6S1/2 state in cm -1 . For s-states the energies were calculated with 5 = —0.20 a.u., for <i-states - 
with 5 = 0.0 a.u., and for p with 5 = —0.09 a.u.. 

DHF DHF+MBPT Experiment [22] 



1 ATlflff 


r^val 


A 






A 
LA 




A 
_A 


65i/2 


0.127368 






0.143085 




0.143099 a 




6pi/ 2 


0.085616 


9163.6 




0.092173 


11172.2 


0.092167 


11178.2 


6^3/2 


0.U83 lob 


9565.3 




0.089609 


1 1 TO A n 

11 1 34.9 


0.U89642 


I 1 TOO yi 

II ^ oz.4 


5d 3 /2 


0.064419 


13815.7 




0.076995 


14503.3 


0.077035 


14499.5 




0.064529 


13791.5 




0.076459 


14621.0 


0.076590 


14597.1 


< s l/2 


U.UOOlo 1 


15841.8 




0.058475 


loODo.U 


U.UOoD^iO 


loOoO.O 


7P1/2 


0.042021 


18731.4 




0.043868 


21773.9 


0.043928 


21765.7 


7P3/2 


0.041368 


18874.8 




0.043041 


21955.4 


0.043103 


21946.7 




"For the . 


ground state E va i = IP (Cs) 


= 31406.71 ■ 


cm -1 33]. 






Rb : 


1 - 20s, 


2 - 


20p, 3 - 19d, 


4 - 19/, 5 


-13(7; 






Cs : 


1 - 23s, 


2 - 


23p, 3 - 23d, 


4 - 26/, 5 


-1%. 





Using these basis sets we solved the multi-particle Shrodinger equation ((7|) and found 
the wave functions of low-lying states. As discussed in [3l[ and demonstrated in [3(1 a 
proper approximation for the effective Hamiltonian can substantially improve an agreement 
between calculated and experimental spectra of multielectron atom. One can introduce an 
energy shift 5 and replace S(£?) — > H(E — 5) in the effective Hamiltonian, Eq. (JHJ). We have 
determined 6 from a fit of theoretical energy levels to experimental spectrum. Using only one 
fitting parameter 8 we reproduced the experimental energies for 12 low-lying states for Li 
and for 10 low-lying states for Na and K with accuracy 0.1-0.2%. To reproduce the low-lying 
energy levels with the same 0.1-0.2% accuracy for heavier Rb and Cs we used three fitting 
parameters (different shifts 5 for different partial waves). An illustrative comparison for the 
heaviest atom Cs (55 electrons) is presented in Table |1J It is worth noting that an empirical 
introduction of shifts 5 mimics higher-order many-body corrections in perturbation theory. 
We will estimate theoretical uncertainty based on sensitivity of our results to variation in 
these shifts. 

IV. QUADRUPOLE AND OCTUPOLE POLARIZABILITIES 

To reiterate major steps of the formalism described in Section m we determined ground 
state wave functions from the effective many-body Shrodinger equation ((Zj), calculated 
dressed electric multipole operators T e s, solved inhomogeneous equation © and computed 
valence parts a\ of dynamic polarizability with Eq. (jlO)) . Additional contributions af of 
core-exited states were calculated using RRPA method. 

Calculation of dynamic polarizabilities with u = gives us the static polarizabilities. We 
provide these data in Tables |H] and IHII and compare them with other results. To estimate 
uncertainties we present in the Tables results of pure DHF calculations and compare them 
with DHF+MBPT ones. The uncertainties of calculations are associated with higher orders 
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TABLE II: Static quadrupole polarizabilities a% for ground states of alkali-metal atoms in a.u. We 
present valence contributions for the cases of pure DHF and DHF+MBPT, and core contributions. 
Final values were determined as sum of a 2 (DHF+MBPT) and a 2 . 





Li 


Na 


K 


Rb 


Cs 


"2 (DHF) 


1485.5 


2230.3 


7049 


9790 


16613 


a v 2 (DHF+MBPT) 


1424.5 


1883.6 


4983 


6488 


10388 


a c 2 (RRPA) 


0.1 


1.5 


16 


35 


86 


Final 


1424(4) 


1885(26) 


5000(45) 


6520(80) 


10470(390) 




Other works 








Patil and Tang [35] 


1393 


1796 


4703 


6068 


10260 


Patil and Tang [14] 


1403 


1807 


4760 


6163 


10400 


Yan et al. [34] 


1423.266(5) 










Marinescu et al. f!3l 


1424 


1878 


5000 


6495 


10462 


Spelsberg et al. [36] 


1423 


1879 


5001 






Maeder and Kutzelnigg 10 


1383 


1799 


4597 


5979 


9478 



of the MBPT which are taken into account only partially. The heavier the atom, the larger 
MBPT contribution is and we expect theoretical accuracy to become worse. For instance, 
the MBPT correction to the static quadrupole polarizability a 2 for Li is only 4%, while for 
Cs it attains 38%. For static octupole polarizabilities the MBPT corrections are larger 
and range from 5% for Li to 48% for Cs. 

Let us turn to estimates of theoretical uncertainty of quadrupole polarizabilities. Es- 
sentially it is based on sensitivity of our results to semiemprically introduced shifts 5. As 
mentioned in Section ILTT1 an introduction of these shifts mimics omitted higher-orders of per- 
turbation theory. We estimate the theoretical error bar as a half of the difference between 
ab initio (5 = 0) value and result with semiempirically chosen 5. Further an overwhelm- 
ing contribution to static 2 z -pole polarizabilities Eq. (j^J) comes from the lowest-lying valence 
state of proper angular symmetry. Since we recover experimental energies almost exactly 
(see Table H}, the theoretical uncertainty is determined by an accuracy of calculation for 
electric-multipole operators of principal transitions. We write 

5a 2 {0) 



(ns\ 


T 2 


n'd) 


& — 


(ns\ 


T 2 \ 


n'd) s=0 


(ns\ 


T 2 \ 
J o 


n'd) 


<5=0 





a 2 (0) 

where ns denotes the ground state and n'd stands for lowest-lying valence d-states. For 
example, following this procedure we obtain an error bar of 0.3% for Li. Our result of 
1424(4) for Li is in excellent agreement with the value 1423.266(5) from benchmark high- 
accuracy variational non-relativistic calculations by Yan et al. |3j|. We estimate theoretical 
uncertainties for octupole polarizabilities to be at 10% level for heavy atoms. Our results 
for static polarizabilities are listed in Tables ITT1 and Hill In these Tables we also compare our 
results with the predictions by other authors. We find that for light atoms there is a good 
agreement between different results except the values obtained by Maeder and Kutzelnigg 
jlDI| are consistently smaller. As the number of atomic electrons increases, the correlation 
effects become more pronounced and discrepancies between results from different groups 
grow larger. Marinescu et al. [ljjj used a model potential with five adjustment parameters 
obtained by fitting to experimental energy levels. Core-polarization was included in the 
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TABLE III: Static octupole polarizabilities 03 for ground states of alkali-metal atoms in 10 4 a.u. We 
present valence contributions for the cases of pure DHF and DHF+MBPT, and core contributions. 
Final values were determined as sum of Q3 (DHF+MBPT) and a|. 





Li 


Na 


K 


Rb 


Cs 


aS (DHF) 


4.185 


6.888 


28.10 


41.50 


76.49 


aS (DHF+MBPT) 


3.957 


5.536 


17.73 


23.66 


39.43 


a§ (RRPA) 





0.001 


0.01 


0.03 


0.10 


Final 


3.957 


5.54 


17.7 


23.7 


39.5 




Other works 










Patil and Tang [35] 


3.871 


5.287 


16.07 


20.73 


33.12 


Patil and Tang [14] 


3.986 


5.430 


16.30 


20.97 


33.33 


Yan et al. [£4] 


3.965049(8) 










Marinescu et al. f!3l 


3.969 


5.552 


17.69 


23.69 


39.53 


Spelsberg et al. [36] 


3.927 


5.486 


19.14 






Maeder and Kutzelnigg [10] 


3.680 


5.117 


15.02 


21.27 


33.99 



pseudo-potential and they also included effects of shielding (or field dressing) in the multipole 
operators. Patil and Tang w|| also used effective potential in their calculations to obtain 
the wave functions of excited states, but they used one-parametric potential and did not 
shielding in the multipole operators. Generally, our results are in a good agreement with 
all results except for values by Maeder and Kutzelnigg . One of possible reasons for this 
discrepancy is that these authors used very small number of basis functions (e.g. only 5 basis 
orbitals for p, d, and / partial waves) while «2 and as polarizabilities are very sensitive to 
details of construction and saturation of basis sets. 

Also shown in Tables |H] and IIHI are the corrections af due to core-excited states. These 
quantities are essentially polarizabilities of singly- charged ions of alkali-metal atoms. Only 
disregarding distortion of the core by the valence electrons, one may identify corrections 
af as core polarizabilities. For static quadrupole polarizabilities their relative contribution 
to the total polarizabilities ranges from 0.01% for Li to 0.8% for Cs. The core corrections 
to static octupole polarizabilities are even smaller (just 0.25% for Cs). Relative smallness 
of af terms for static polarizabilities may lead one to a wrong assumption that the core 
excitations may be disregarded in calculations of van der Waals coefficients C n . In fact the 
expression (JHJ) for C n contains integration over an infinite range of frequencies u. While the 
region around to = does provide the dominant contribution to C n , the high-frequency tail 
of the polarizability is still important. As uj —>■ 00 the core polarizability overpowers valence 
contribution. In fact, one of the points of the paper p] was to explicitly demonstrate that 
for heavy atoms the core polarizability may contribute as much as 15% to Cq dispersion 
coefficient. Here using RRPA calculations of af(iu) core polarizability we will arrive at a 
similar conclusion for higher-multipole coefficients Cs and C\q. 

We calculated the core polarizabilities in the framework of relativistic random-phase ap- 
proximation method (RRPA). Essentially we extended approach of Johnson et al. [2^ and 
incorporated frequency dependence into the calculations. Compared to Ref. [26] we also em- 
ployed a different numerical technique using B-spline basis sets. With our newly developed 
code we recover the previous results for static dipole and quadrupole polarizabilities. We 
found that unusually large basis sets of 100 B-splines were required to achieve a numerical 
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TABLE IV: Static quadrupole polarizabilities 0^(0) of singly-charged ions of alkali-metal atoms 
(core polarizabilities). Results marked RRPA are results of our calculations; these numerical values 
are identical to those by Johnson et al. [2H]. All values are in atomic units. 

Li+ Na+ K+ Rb+ Cs+ 

RRPA Oil L52 io 351 861 

Patil [2L .ag] 1.64(15) 18.2(3.0) 42(3) 128(40) 

Freeman and Kleppner [£9j 1.91(15) 



TABLE V: Static octupole 


polarizabilities 


a§(0) of sinj 


dy-charged ions 


of alkali-metal atoms (core 


polarizabilities). 


All values 


are in atomic units. 








Li+ 


Na+ 


K+ 


Rb+ 


Cs+ 


This work 


0.17 


7.5 


110 


314 


1014 


Patil [22] 






95(10) 


280(40) 


1220(200) 



convergence, especially for octupole polarizabilities of heavy atoms. Finally, we present a 
comparison of the computed RRPA static quadrupole and octupole core polarizabilities with 
other works in Tables IIVI and Patil |37l . l38| has inferred these polarizabilities analyzing 
Rydberg spectra of alkalis. His results are in a uniform agreement with our ab initio values. 



V. VAN DER WAALS COEFFICIENTS 

From general formula (JHJ) dispersion coefficients may be expressed as 

Cft = C a b(l, 1), 

C? = C a6 (l,2)+C o6 (2,l), 

CT fe = C ab (2, 2) + C ab (l, 3) + C ab (3, 1) . (11) 

Here the coefficients C ab (l,l') are quadratures of atomic 2 l — and I 1 ' — pole dynamic polariz- 
abilities 

3 /-oo 

(12) 
(13) 
(14) 
(15) 

Calculations of dynamic polarizabilities were discussed in the previous section and here we 
proceed to evaluation of the dispersion coefficients. 

The computed C$ and C\o coefficients for homonuclear and heteronuclear species are 
presented in Tables IVIT IIXI The dispersion coefficients Cq were tabulated previously in 



C a6 (i,i) 


3 

7T . 


POO 

la 




iuS) 


«?( 


iu)du, 


c; 6 (i,2) 


15 
~ 2^ 


POD 

Jo 


a\ 


(iu] 


)a b 2 


(iu)du. 


C a fe(2, 2) 


35 

7T 


rocs 

Jo 


a a 2 


{iu) 




[iu)du. 


C a6 (l,3) 


14 

71 


POO 

Jo 


al 


(iu] 


)«3 


(iu)du 
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Refs. |6|, [17(. This completes the first application of relativistic many-body methods of 
atomic structure to calculations of leading long-range interactions between ground-state 
alkali-metal atoms. 

To estimate uncertainties in our values we notice that the main value of the quadratures, 
Eqs. f|12jl -(|15 jl is accumulated in the low-frequency region uj « 0. Therefore the error may 
be expressed via uncertainties in the static multipole polarizabilities 



SC ab (l,l') = I /M0)V ( Sa v {0) 



1/2 



induced in Cg b is 



The required uncertainties Sai (0) were estimated in Section |IV] and Ref. p. The error 

1/2 



SCf = {(SC ab (l, 2)) 2 + (5C ab (2, l)) 2 } 



Here we assumed that a ^ b. The formulas for homonuclear dimers may be derived in a 
similar manner. The resulting theoretical uncertainties for C$ coefficients range from 0.5% 
for Li 2 to 4% for Cs dimer. We anticipate uncertainty in Cio coefficients to be better than 
10%. 

It is instructive to consider the effect of core excitation contribution af(iu;) to dynamic po- 
larizabilities and thus to C„ coefficients. Such corrections are omitted in the model potential 
calculations such as Ref. [3, El. To illuminate the relative contributions of core-excitations 



we computed C n coefficients by keeping only the valence contributions to the total dynamic 
polarizabilities 

aii(iu) — > a v t {iuj). 



Such calculated dispersion coefficients are marked as Cg and in Tables IVIHIXj while 
values marked "final" were obtained with an additional inclusion of core excitations. Com- 
paring these values, we observe that relative contribution of af(iu) term grows rapidly as 
the number of atomic electrons increases. For example, examining Table IVD we see that core 
correction to Cg for Li is only 0.2%, while for Cs it is 10%. For Cio coefficients the core 
contributions for all atoms are slightly smaller. Still for Cs core excitations contribute 8% 
to the Cio coefficient. 

A comparison with results by other authors is presented in Tables fVTHlXI There is good 
agreement for light Li and Na atoms. For heavier atoms, in particular for Cs, there is 
discrepancy at the level of 10% for Cg and 20% for Cio coefficients. Such tendency may be 
attributed to two factors. First, correlations become enhanced for heavier atoms. Another 
cause is that model-potential calculations such as Ref. [HQ disregard contribution of core- 



excited states. This corresponds to the valence term denoted as C^ in Tables IVIHIXI As 
mentioned above the core-excited states contribute at the level of 10% for Cs. If we disregard 
this contribution, we see that the model-potential results are in a reasonable agreement with 
our C^ values. 

Only very recently interpretation of experiments with ultracold atoms allowed several 
groups to reduce uncertainties in the C$ coefficients to a fraction of a per cent (H, 0, H|. 
These inferred coefficients are in an excellent agreement with our values predicted using 
many-body perturbation theory |6| . Even more refined understanding of details of ultracold 
collisions led very recently to constraints on higher-multipole coefficient Cg for Rb 2 0, HI 
and Cs dimer . In Table IVII we present a comparision with these inferred values. Our 
computed value for Rb 2 5.77(8) agrees well with C 8 = 5.79(49) by van Kempen et al. Q and 
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TABLE VI: van der Waals Cg coefficients in 10 5 a.u. for homonuclear dimers. Cg values include 
only valence contributions. The final values were determined as combination of DHF+MBPT 
method for valence contributions with RRPA calculations for core excitations. 





Li 


Na 


K 


Rb 


Cs 


u 8 


0.832 


1.15 


a nn 


O.O 1 


n ifi 
y . iu 


Hi n q 1 


0.834(4) 


1.160(18) 


4 90('5 N I 


5 77f8"l 


1 9(A) 




Other theoretical works 








Patil and Tans' [141 

A CLVLl till LI iOjllci | J-t: 


0.8183 


1.090 


3.892 


5.258 


9.546 


Yan et al. [£4] 


0.834258(4) 










TVTan'nPKPn pf al 13 

J.V±CXL lllCDb LL C lj (XI . I 1U 


0.8324 


1.119 


4.096 


5.506 


9.630 


*■» t~i o I o Kovrf nf q -ir\ 


0.8303 


1.141 


A 01 1 






A/I fi pn pr ann T\ n f'zpl n l cecr 10 

±V±OjC'U.C'1 J. V Li IjZjCllllggj |XVJ 


0.8089 


1.098 


3 834 


5.244 


9 025 




Experiment 








van Kempen et al. [7J 








5.79(49) 
6.09(7) 




Marte et al. [8] 








5.73 




Leo et al. [9] 










8.4(4) 



Cg = 5.73 by Marte et al. [3]. However, we disagree with 1%-accurate result 0] of 6.09(7) 
by four standard deviations. This 1%-accurate result was obtained in Ref. [2J by setting 
additional constraints on the singlet potential of Rb dimer while including higher-multipole 
van der Waals coefficients C\\ and Cyi in the fit. For Cs2 the inferred value by Leo et al. 
is C 8 = 8.4(4), it disagrees with our prediction, 10.2(4) by more than four standard 
deviations. It is worth noting that while for Rb the inferred value lies above our result, for 
Cs the situation is reversed and our value is larger. 

To conclude, we calculated static and dynamic quadrupole and octupole polarizabilities 
for ground states of Li, Na, K, Rb, and Cs atoms. The calculations were carried out 
using accurate relativistic many-body methods of atomic structure. With the computed 
polarizabilities we evaluated Cg and C\q van der Waals coefficients for homonuclear and 
heteronuclear dimers and estimated theoretical uncertainties. The estimated uncertainties 
for Cg coefficients range from 0.5% for Li 2 to 4% for Cs2- We have highlighted the role of 
usually omitted core excitations in calculation of Cg, and Cio coefficients and found that 
their contribution is important for heavy atoms K, Rb, and Cs. 

We would to thank Mikhail Kozlov for comments on the manuscript. This work was 
supported in part by the National Science Foundation. The work of S.G.P. was additionally 
supported by the Russian Foundation for Basic Research under Grant No. 02-02-16837-a. 



[1] A. Dalgarno and W. D. Davison, in Adv. At. Mol. Phys., edited by D. Bates and I. Estermann 

(Academic Press, New York, 1966), vol. 2, pp. 1-32. 
[2] J. Weiner, V. S. Bagnato, S. Zilio, and P. S. Julienne, Rev. Mod. Phys. 71, 1 (1999). 
[3] J. L. Roberts, J. P. Burke, Jr., N. R. Claussen, S. L. Cornish, E. A. Donley, and C. E. Wieman, 

Phys. Rev. A 59, 578 (1999). 
[4] P. J. Leo, C. J. Williams, and P. S. Julienne, Phys. Rev. Lett. 85, 2721 (2000). 



10 



TABLE VII: van der Waals C\q coefficients in 10 a.u. for homonuclear dimers. Cf values include 
only valence contributions. 





Li 


AT 

JNa 


K 


Kb 


Cs 


°10 


0.734 


1.12 


5.18 


7.55 


14.7 


b mal 


0.735 


1.13 


5.37 


7.96 


15.9 




Ut her 


works 








Patil and Tang [14] 


0.7289 


1.068 


4.789 


6.833 


13.58 


Yan et al. [34] 


0.73721(1) 










Marinescu et al. [13J 


0.7365 


1.107 


5.248 


7.665 


15.20 


Spelsberg et al. [26J 


0.7306 


1.113 


5.431 






Maeder and Kutzelnigg [10J 


0.6901 


1.036 


4.522 


6.836 


13.01 



TABLE VIII: van der Waals Cs coefficients in 10 5 a.u. for heteronuclear dimers. C| values include 
only valence contributions. 





Li-Na 


Li-K Li-Rb 


Li-Cs Na-K Na-Rb 


Na-Cs 


K-Rb K-Cs 


Rb-Cs 


Ci 


0.982 


1.91 2.26 


3.07 2.18 2.56 


3.43 


4.64 6.13 


7.04 


Final 0.988(1L 


) 1.95(2) 2.34(3) 


3.21(10) 2.24(3) 2.66(4) 


3.62(12) 


4.93(6) 6.62(19) 


7.69(22) 








Other theoretical works 








M 


0.949 


1.852 2.190 


3.049 2.082 2.444 


3.355 


4.531 6.162 


7.111 


[13J 


1.068 


2.517 3.137 


4.586 2.614 3.250 


4.727 


5.123 7.547 


8.120 


[36] 


0.978 


1.911 


2.174 








[5] 


C. Amiot, 


O. Dulieu, R. F. 


Gutterres, and F. Masnou-Seeuws, 


Phys. Rev. A 66, 052506 




(2002). 












[6] 


A. Derevianko, W. R. Johnson, M. S. Safronova, and J. 


F. Babb, 


Phys. Rev. Lett. 


82, 3589 



(1999). 

[7] E. G. M. van Kempen, S. J. J. M. F. Kokkelmans, D. J. Heinzen, and B. J. Verhaar, Phys. 

Rev. Lett. 88, 093201 (2002). 
[8] A. Marte, T. Volz, J. Schuster, S. Diirr, G. Rempe, E. G. M. van Kempen, and B. J. Verhaar, 

Phys. Rev. Lett. 89, 283202 (2002). 



TABLE IX: van der Waals C\q coefficients in 10 a.u. for heteronuclear dimers. Cf values include 
only valence contributions. 





Li-Na 


Li-K 


Li-Rb 


Li-Cs Na-K Na-Rb 


Na-Cs 


K-Rb 


K-Cs 


Rb-Cs 


°10 


0.912 


2.07 


2.55 


3.73 2.48 3.04 


4.40 


6.3 


8.9 


10.6 


Final 


0.916 


2.10 


2.61 


3.84 2.53 3.13 


4.55 


6.6 


9.4 


11.3 










Other theoretical works 










[14] 


0.8859 


1.949 


2.356 


3.379 2.303 2.773 


3.948 


5.724 


8.077 


9.629 


[13J 


0.982 


2.651 


3.413 


5.303 2.949 3.784 


5.844 


6.726 


10.37 


11.79 


[36] 


0.9058 


2.139 




2.553 











11 



[9] P. Leo, C. J. Williams, E. Tiesinga, and P. S. Julienne (2003), private communication. 
[10] F. Maeder and W. Kutzelnigg, Chem. Phys. 42, 95 (1979). 
[11] J. M. Standard and P. R. Certain, J. Chem. Phys. 83, 3002 (1985). 
[12] B. Bussery and M. Aubert-Frecon, J. Chem. Phys. 82, 3224 (1985). 
[13] M. Marinescu, H. R. Sadeghpour, and A. Dalgarno, Phys. Rev. A 49, 982 (1994). 
[14] S. H. Patil and K. T. Tang, J. Chem. Phys. 106, 2298 (1997). 

[15] S. H. Patil and K. T. Tang, Asymptotic Methods in Quantum Mechanics: Application to 
Atoms, Molecules and Nuclei (Springer- Verlag, Berlin Heidelberg, 2000). 

[16] D. A. Varshalovich, A. N. Moskalev, and V. K. Khersonsky, Quantum Theory of Angular 
Momentum (World Scientific, Singapore, 1988). 

[17] A. Derevianko, J. F. Babb, and A. Dalgarno, Phys. Rev. A 63, 052704 (2001). 

[18] A. Derevianko and S. G. Porsev, Phys. Rev. A 65, 053403 (2002). 

[19] V. A. Dzuba, V. V. Flambaum, and M. G. Kozlov, Pis'ma Zh. Eksp. Teor. Fiz. 63, 844 (1996), 

[JETP Lett. 63 882, (1996)]. 
[20] V. A. Dzuba, V. V. Flambaum, and M. G. Kozlov, Phys. Rev. A 54, 3948 (1996). 
[21] V. A. Dzuba, M. G. Kozlov, S. G. Porsev, and V. V. Flambaum, Zh. Eksp. Teor. Fiz. 114, 

1636 (1998), [Sov. Phys.-JETP 87 885, (1998)]. 
[22] S. G. Porsev, Y. G. Rakhlina, and M. G. Kozlov, J. Phys. B 32, 1113 (1999). 
[23] S. G. Porsev, Y. G. Rakhlina, and M. G. Kozlov, Phys. Rev. A 60, 2781 (1999). 
[24] M. G. Kozlov and S. G. Porsev, Eur. Phys. J. D 5, 59 (1999). 
[25] D. Kolb, W. R. Johnson, and P. Shorer, Phys. Rev. A 26, 19 (1982). 
[26] W. R. Johnson, D. Kolb, and K. N. Huang, At. Data Nucl. Data Tables 28, 333 (1983). 
[27] R. M. Sternheimer, Phys. Rev. 80, 102 (1950). 
[28] A. Dalgarno and J. T. Lewis, Proc. Roy. Soc. 233, 70 (1955). 

[29] V. F. Bratsev, G. B. Deyneka, and I. I. Tupitsyn, Bull. Acad. Sci. USSR, Phys. Ser. 41, 173 
(1977). 

[30] S. G. Porsev, M. G. Kozlov, Y. G. Rakhlina, and A. Derevianko, Phys. Rev. A 64, 012508 
(2001). 

[31] M. G. Kozlov and S. G. Porsev, Opt. Spectrosk. 87, 384 (1999), [Opt. Spectrosc. 87 352, 
(1999)]. 

[32] S. G. Porsev, M. G. Kozlov, and Y. G. Rakhlina, Pis'ma Zh. Eksp. Teor. Fiz. 72, 862 (2000), 

[JETP Lett. 72 595, (2000)]. 
[33] C. E. Moore, Atomic energy levels, vol. I-III (NBS, National Standards Reference Data Series 

- 35, U.S. GPO, Washington, D.C., 1971). 
[34] Z.-C. Yan, J. F. Babb, A. Dalgarno, and G. W. F. Drake, Phys. Rev. A 54, 2824 (1996). 
[35] S. H. Patil and K. T. Tang, Chem. Phys. Lett. 301, 64 (1999). 
[36] D. Spelsberg, T. Lorenz, and W. Meyer, J. Chem. Phys. 99, 7845 (1993). 
[37] S. H. Patil, J. Phys. B 27, 601 (1994). 
[38] S. H. Patil, Phys. Phys. A 33, 90 (1986). 

[39] R. R. Freeman and D. Kleppner, Phys. Rev. A. 14, 1614 (1976). 



12 



